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We explore several microscopic mechanisms for breaking the n — fourfold Landau level de- 
generacy in a single-layer graphene. Valley-scattering random potential, Zeeman interaction, and 
electron-phonon coupling are considered in the presence of SU(4)-symmetric Coulomb exchange 
interaction. Among all the mechanisms considered, it is the electron-phonon coupling combined 
with the Zeeman interaction which leads to the full splitting of the n — Landau levels. A recent 
controversy of "valley-first" or "spin-first" breaking of SU(4) symmetry of the n = graphene Lan- 
dau level is examined in light of our results. Existence of midgap states between Landau levels of 
opposite valley polarity are demonstrated. 
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Spin and valley degeneracy of a single-layer graphene 
sheet subject to a perpendicular magnetic field gives rise 
to a quantized Hall conductance [l|, 0] o'xy that changes 
in multiples of four units of conductance quantum /h: 
O'xy = 4:{e'^/h){n + 1/2), n=integer. The conductance 
changes from (taking e'^ /h = 1) —2 to +2 as the chem- 
ical potential passes through the central, n = Lan- 
dau level (LL) which is fourfold degenerate. When the 
strength of the magnetic field increases, these fourfold de- 
generate LL's split in energy into two sublevelsjl, 0, @|j 
and accordingly the Hall conductance steps occur at 
—2,0, and +2. At even higher fields the fourfold de- 
generacy is broken completely, leading to the sequence 
a,j, = -2, -1,0,1,2:3, 4, 5]. 

The manner of the breaking of fourfold degeneracy 
of the central LL has been discussed in a number of 
papers[6, 7, i, i, M, Ell, [13, E, [II- The SU(4) sym- 
metry of the pristine ?i = graphene LL is broken down 
either spontaneously from interaction effects, or due to 
various weak symmetry-breaking terms such as lattice ef- 
fects, Zeeman splitting, etc. Since the degeneracy arises 
from SU(2) symmetry of the electron spin and another 
SU(2) symmetry of the valley, the main focus of discus- 
sion has been whether the spin-symmetry or the valley- 
symmetry breaking occurs first. Abanin, Levitov, and 
Lee argued that the spin symmetry breaking should oc- 
cur first, turning the edge of the Hall sample into con- 
ducting channels0. An experiment carried out shortly 
thereafter seems to confirm their picture j4[. On the other 
hand, subsequent transport experiments [la| which found 
divergent longitudinal resistance seem to rule out the ex- 
istence of such gapless edge states, and the issue of "spin- 
first" or "valley-first" symmetry breaking appears by no 
means settled. 

In this paper, we revisit the "hierarchy problem" of 
the central LL splitting in a single-layer graphene within 
the self-consistent Hartree-Fock theory, while considering 
several SU(4)-symmetry breaking terms explicitly. Fol- 
lowing the general approach, we adopt the continuum 
description of the graphene dynamics using the spinor 



Spin (cr =t,i) and valley index 



T = ± are introduced to classify the spinors formed from 
a- and 5-sublattice electrons. The Landau level problem 
with the perpendicular magnetic field can be treated by 
the Hamiltonian 
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With the non-commuting operators obeying — 
ifi^ /Pg, Ib = ^Ti/eB being the magnetic length, one can 
form a set of canonical operators a = [Ib / ^/2K){px + 



ipy), a+ = {lB/V2h){px - ipy), [a, a+] = 1. A cyclotron 
frequency uj — ^/2vf/Ib {vF=Feiim velocity) has been 
introduced above. 

In writing down the Hamiltonian in the manifestly 
SU(4)-symmetric form above, we have implemented the 
rotation of the r = — spinor, ipa-^ — > cryipa--- Using the 
complete set of normalized eigenfunctions given by 
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one may expand the field operator as ipariT) — 
Y.n,mXnmT{r)"fnraar- Here 0„m is the oscillator wave 
function, 771 is the guiding center coordinates, and Xnmr 
equals Xnm defined in Eq. ([2]) if r = +, but equals ayXnm 
when r = — . 

As the primary interest of this paper is in understand- 
ing the mechanism of level splitting within the central LL, 
we carry out the projection to 7i = LL states. The ki- 
netic energy gets completely quenched, whereas Coulomb 
interaction within this LL reads 
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The reference to the LL index n = has been dropped. 
The summation over the repeated guiding center coordi- 
nates is implicit. For numerical purpose, we work with a 
torus geometry of dimension L^; x Ly and use the Landau 
gauge for which the wave functions are 



4>ni{r) = 



1 



^'y--e~"^^y-y^^^\yra='^m. (4) 
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The Coulomb Hamiltonian in this basis reads 
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where the Fourier-transformed Coulomb potential 
V{k) = / (PrV{r)e'^''''^ is shown. Whereas kx,ky run 
over all integer multiples of 2'it/Lx and 27r/Ly, the guid- 
ing center coordinates mi and m2 span 1 through N^, 
the number of flux through the lattice given by 
— LxLy. 

Having established a discretized Hamiltonian, we 
solve it within the Hartree-Fock theory using the self- 
consistent parameter 



A 



(mi,m2) = (7. 
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with an arbitrary pair of guiding center indices mi, 7712 
and the spin-valley indices. For the reason that Hartree 
term offers only a chemical potential shift for the uni- 
form solutions we find, and that Hartree interaction does 
not break the SU(4) symmetry, we will be exclusively 
concerned with the exchange Hamiltonian, H^^. 

Among the possible SU(4) symmetry-breaking terms 
we consider the following three: (i) Zeeman field: 
= Bcr J2maT '^^niarlniar, (u) Valley-scattcring im- 
purity: iJ™P = Z)™^^ KiT7+^^7m<Tr- We take as 
a random number of width W: Vm & [—W/2,W/2], 
and (iii) Valley-scattering electron-phonon coupling[16]: 

—t). Projected onto the central LL and treated in the 
mean-field manner, this last Hamiltonian becomes 



mar 



where U' = ?7/v27r[l7|. For convenience, we use U in- 
stead of U' from now. 

The total Hamiltonian we will consider is H = + 



cl — ph 
MF ■ 



Although all these terms indi- 



vidually may have been analyzed in various ways in the 
past, it is our belief that their combined effects and possi- 
ble competition among different symmetry-breaking ten- 
dencies have never been studied in the presence of the 



FIG. 1: (color online) Plot of energy levels at half-filling 
obtained from self-consistent solutions of H^'^ (black solid 
lines), H^^ + (red dash lines), H^^ + iJ^P (magenta 
crosses), and H^^ + + H'"'"" (blue filled squares). The 
impurity broadening is W — 0.1 and the Zeeman field is 
Bcr ~ 0.08. The system size used is A^^ = 50. 



Coulomb exchange interaction in a self-consistent man- 
ner. We made extensive numerical simulation at zero 
temperature to identify which of the combinations of the 
Zeeman, impurity, and electron-phonon coupling-induced 
interactions would lead to the full splitting of the four- 
fold degeneracy. Both quarter-filled and half-filled cases 
were examined. The 3/4-filled case can be deduced by 
symmetry from the results of I/4-filled case. 

With H^-^ alone and at half-filling, the initial four- 
fold degeneracy of the LL is split into two sublevels with 
energies at ±E^'^, where the scale E^-^ is set by the 
exchange energy interaction. In our convention, c^/kI 
(K=dielectric constant) is taken to unity, and in such a 
unit we obtain E^-^ w 0.5. For the quarter-filled case, a 
similar situation arises with one LL at an energy below 
the chemical potential and three degenerate LL's whose 
energy lies above it. From these exercises we learn that 
the full energy splitting of the central LL requires more 
than the Coulomb exchange effect alone. Which of the 
spin and valley symmetry remains intact is completely 
arbitrary at this point. 

Still at half-filling, inclusion of the Zeeman field to 
the Coulomb exchange now ensures that the symmetry 
breaking occurs along the spin direction, with the LL 
energies at ±{E^^ + B<t)- The valley-SU(2) symmetry 
is preserved under the addition of H^. It thus appears 
that more than one SU(4) symmetry breaking mecha- 
nism need to be present to fully split the degeneracy. We 
find that further inclusion of the valley-scattering impu- 
rity, H'^^ + + ff™P, does not result in additional 
splitting of the levels. The previously twofold degenerate 
states for each guiding center m undergoes splitting by 
±Vm, and give rise to broadened energy levels of width 
W. The numerically obtained energy levels for several 
combinations of terms at half-filling can be found in Fig. 

m 

Actually, the Coulomb exchange Hamiltonian per- 
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turbed by two kinds of Zeeman fields separately act- 
ing on the spin and the valley spaces, as in H^-^ + 

BaT,rnar'^l^arlrnar + Br T,mar '^Imar^/mcrr , doCS ex- 
hibit a full lifting of the fourfold degeneracy with the en- 
ergies given at (£'^^-hBM)±Sm and at -{E^^+BM)±Bm. 
Here Bm and Bm refer to the larger and the smaller of 
the two Zeeman fields, respectively. One can also see 
that the electron-phonon coupling iJ^p''^ provides the 
required valley Zeeman field, acting along the x-axis of 
the valley spin. 
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FIG. 2: (color online) (a) Phase diagram of the central LL 
at half-filling with varying Zeeman field {2Ba) and electron- 
phonon coupling strength (U). In region I, only the spin de- 
generacy is broken. In region II, both spin and valley symme- 
tries are lost, but the main level splitting takes place in the 
valley direction. A first-order phase boundary separates the 
two regions, (b) Schematic energy levels in each region. The 
energy level separations and spin- valley quantum numbers for 
each LL are specified. 

Aided by these ideas, we next consider + + 



^MF^ ■ Figure [2] shows the phase diagram for such a 
model, at half-filling, spanned by two interaction param- 
eters {U,2Ba)- There are two phases found here, called 
I and II, distinguished by the number of levels split. In 
region I, where the Zeeman effect dominates, only one, 
spin-polarized level splitting is observed. The valley- 
splitting order parameter Ao--,ct+(to, m) becomes zero in 
this region, still preserving the valley symmetry. In the 
[/-dominated region II, the main polarization direction 
is along the valley axis, and the Zeeman field contributes 
to the sublevel splitting equal to 2Bc,. Here indeed, the 
full breaking of the fourfold degeneracy is obtained. The 
phase boundary taking place along 2B^ w [/ in the phase 
diagram is first-order. 

The phase diagram for quarter-filling is similar, as 
shown in Fig. [3) The phase boundary now taking place 
exactly at 2Bcr = U separates the "spin- first" -split region 
I from the "valley-first" -split region II. The two centrally 
located LL's cross in energy at the phase boundary. For 
instance, ^^(t, — ) < E{i, +) energy hierarchy in region I 
crosses over to E{^,—) > E{1,+) in region II. We have 
checked that the inclusion of the impurity does not alter 
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FIG. 3: (color online) Phase diagram of the central LL at 
quarter-filling. Meaning of the symbols are the same as in 
Fig. [2] Only one LL lies below the Fermi level here. 



the basic features of the phase diagram shown in Figs. 
and [3] as long as W remains small compared to both U 
and 2B„. 

The relevant energy scales 2B^ and [/ in a graphene 
layer are comparable, as recently discussed in Ref. 
[16l |. Both energy scales are of order Ba in units of 
[Kelvin/Tesla], and therefore it should be quite possible 
that graphene samples with either 2Ba > U or 25^ < U 
exist. Then according to Ref. 0], the edge of the half- 
filled graphene quantum Hall system can be either con- 
ducting {U < 2B„) or insulating {U > 2B„). Another 
interesting possibility suggested by our search is the tran- 
sition between the two LL splitting scenarios driven by 
the relative strengths of electron-phonon coupling and 
Zeeman energies. We speculate that applying a mechan- 
ical pressure, such as stretching, to the graphene will 
influence U without changing the Zeeman energy, and 
might allow one to probe the phase transition between 
regions I and II. The bond-CDW order associated with 
region II at half-filling and regions I and II for quarter- 
filling should leave a mark in the electronic spectrum, 
which can be probed by STM. 

A Landau level with a quantum number (cr, r) occurs at 
the same energy as another state with (cr, r), whereas the 
same is not true with (ct, r) due to the Zeeman splitting. 
One can then imagine domain walls separating the two 
LL states with opposite valley polarities in a macroscopic 
sample. Assuming the spinless case for simplicity, the 
physics of such a domain wall can be captured in a set of 
differential equations 

i{dy+yk)[u-+v+]+m^{y)[u+-v^] = -ek[u+-v^] 

i(dy-yk)[u+-vJ\-m^{y)[u^+v+\ = -ek[ii-+v+\ 

i{dy+yk)[u--v+]+m^{y)[u+ + vJ\ = ek[u++vJ\ 

i{dy-yk)[u++vJ\-m^{y)[u^-v+\ = efe[u_-w+]. (8) 

We have written the eigenfunction associated with the a- 
and 6-sublattice as {u■T,v^■)^ also distinguished by their 
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FIG. 4: (color online) The energy Sk obtained from Eq. @ 
with rriv = 0.9. 



valley index r, used the linear gauge {Ax — —By, Ay — 
0) , and taken out the a;-dependence of the wave function 



as {Ur, 



e (ur(y), fx (?/))• The y-dependent mass 



gap due to the inter- valley scattering is written m^ijj), 
and Uk abbreviates y — k. 

As an example of the influence of the sign change of 
the mass gap on the energy spectra, consider the case 
with an abrupt sign change as in m^{y) = mvSgn(y). 
Away from y — Q the mass gap is uniform, and a pair of 
solutions with Sk — is found for w+ — — 0, and 

U- = ±u_|_ = (j){y — k) respectively, where 4i{y — k) is the 
Gaussian function peaked &ty = k. Since k is still a good 
quantum number, the eigen energies Sk can be solved for 
each k separately. By matching the wave functions at 
?/ = we derive the following equation determining the 
energies. 



m^-Ek _ H(pk,k)H'(pk,-k) 
rriy+ek H{pk,-k)H'{pk,k)' 



(9) 



where pk = {e\ — m1)/2 < 0, and H{pk,k) is Hermite 
polynomial of negative order pk ■ Self-consistently solving 
the equation for each k gives rise to the energy band 
shown in Fig. [H 

The level crossing predicted here will be of particu- 
lar relevance if the chemical potential should lie between 
the two valley-split LL's. In such a case, a pair of gap- 
less one-dimensional channels will cross the Fermi level, 
similar to the edge channel in the spin-split case first dis- 
cussed in Ref. Unlike the edge channels, the metallic 
channel predicted here can be formed at the bulk when- 
ever a domain boundary separates the opposite valley 
states. Referring to our phase diagram in Fig. [21 the two 
LL's lying closest to the Fermi level always carry oppo- 
site spins, therefore a domain wall state connecting the 
occupied and the unoccupied LL's will have to involve 
twist in both spin and valley. On the other hand, the 
quarter-filled case (region I) offer a better chance for ob- 
serving the domain wall between two spin-polarized LL's 
which differ only in their valley polarities. We therefore 
suggest that the "spin-first" split LL with axij = — 1 is 
the most likely platform to observe the metallic domain 



walls. 

In summary, we explored the hierarchy issue of the 
central LL symmetry breaking of a graphene layer. Self- 
consistent Hartree-Fock theory was employed, taking into 
consideration several kinds of SU(4) symmetry breaking 
terms. The competitive nature of the valley-splitting 
(due to electron-phonon interaction) and spin-splitting 
(due to Zeeman interaction) leads to a phase diagram 
with either "spin-first" or "valley-first" level splitting. 
Existence of a new kind of gapless state when LLs with 
opposite valley polarities form a domain wall is demon- 
strated. 

H. J. H. is supported by the Korea Science and Engi- 
neering Foundation (KOSEF) grant funded by the Ko- 
rea government (MEST) (No. ROl-2008-000-20586-0), 
and in part by the Asia Pacific Center for Theoretical 
Physics. One of us (H. J. H.) wishes to thank P. Kim for 
hospitality during the completion of this work and for 
discussion. 



[2: 

[3: 

[4; 

[s: 
[6: 
[7: 

[8 
[9 
[W 

[11 

[12 

[13 
[14 

[15 

[16 

[17 



Electronic address: hanjh@skk u.edu] 

K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, 

M. I. Katsnelson, I. V. Grigorieva, S. V. Dubonos, and 

A. A. Firsov, Nature (London) 438, 197 (2005). 

Y. Zhang, Y.-W. Tan, H. L. Stormer, and P. Kim, Nature 

(London) 438, 201 (2005). 

Y. Zhang, Z. Jiang, J. P. SmaU, M. S. Purewal, Y.-W. 
Tan, M. FazloUahi, J. D. Chudow, J. A. Jaszczak, H. L. 
Stormer, and P. Kim, Phys. Rev. Lett. 96, 136806 (2006). 
Dmitry A. Abanin, Kostya S. Novoselov, Uli Zeitler, 
Patrick A. Lee, A. K. Geim, and L. S. Levitov, Phys. 
Rev. Lett. 98, 196806 (2007). 

Z. Jiang, Y. Zhang, H. L. Stormer, and P. Kim, Phys. 
Rev. Lett. 99, 106802 (2007). 

K. Nomura and A. H. MacDonald, Phys. Rev. Lett. 96, 
256602 (2006). 

Dmitry A. Abanin, Patrick A. Lee, and Leonid S. Levi- 
tov, Phys. Rev. Lett. 96, 176803 (2006). 
M. O. Goerbig, R. Moessner, and B. Doucot, Phys. Rev. 
B 74, 161407(R) (2006). 

J. Alicea and M. P. A. Fisher, Phys. Rev. B 74, 075422 

(2006) . 

Igor F. Herbut, Phys. Rev. B 75, 165411 (2007). 

J. N. Fuchs and P. Lederer, Phys. Rev. Lett. 98, 016803 

(2007) . 

L. Sheng, D. N. Sheng, F. D. M. Haldane, and Leon 
Balents, Phys. Rev. Lett. 99, 196802 (2007). 
Motohiko Ezawa, J. Phys. Soc. Jpn. 76, 094701 (2007). 
V. P. Gusynin, V. A. Miransky, S. G. Sharapov, and I. 
A. Shovkovy, Phys. Rev. B 77, 205409 (2008). 
J. G. Checkelsky, Lu Li, and N. P. Ong, Phys. Rev. Lett. 
100, 206801 (2008); Phys. Rev. B 79, 115434 (2009). 
Kentaro Nomura, Shinsei Ryu, and Dung-Hai Lee, 
larXiv:0906.0 159vl. 

A more general valley symmetry-breaking order parame- 
ter AaT,cr-{m,m') can be chosen with m 7^ m'. Perhaps 
due to the localized nature of LL wave functions, actual 
self-consistent values are negligibly small unless m = m' . 



